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ABSTRACT 



Context. Results of forward modelling of acoustic wave propagation in a realistic solar sub-photosphere with two cases of steady 
horizontal flows are presented and analysed by the means of local helioseismology. 

Aims. This paper is devoted to an analysis of the influence of steady flows on the propagation of sound waves through the solar 
interior. 

Methods. The simulations are based on fully compressible ideal hydrodynamical modelling in a Cartesian grid. The initial model 
is characterised by solar density and pressure stratifications taken from the standard Model S and is adjusted in order to suppress 
convective instability. Acoustic waves are excited by a non-harmonic source located below the depth corresponding to the visible 
surface of the Sun. Numerical experiments with coherent horizontal flows of linear and Gaussian dependences of flow speed on depth 
are carried out. These flow fields may mimic horizontal motions of plasma surrounding a sunspot, differential rotation or meridional 
circulation. An inversion of the velocity profiles from the simulated travel time differences is carried out. The inversion is based on 
the ray approximation. The results of inversion are then compared with the original velocity profiles. 

Results. Results of forward modelling of acoustic wave propagation in a realistic solar sub-photosphere with two cases of steady 
horizontal flows are presented. The influence of steady flow on the propagation of sound waves through the solar interior is analysed. 
A time-distance analysis technique is applied to compute the direct observable signatures of the background bulk motions on travel 
times and phase shifts. This approach allows direct comparison with observational data. Further, we propose a method of obtaining 
the travel-time differences for the waves propagating in sub-photospheric solar regions with horizontal flows. The method employs 
directly the difference between travel-time diagrams of waves propagating with and against the background flow. 
Conclusions. The analysis shows that the flow speed profiles obtained from inversion based on the ray approximation differ from the 
original ones. The difference between the original and observed profiles is caused by the fact that the wave packets propagate along 
the ray bundle, which has a finite extent, and thus reach deeper regions of the sub-photosphere in comparison with ray theory. 
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1. Introduction 



Recent observational studies in local time-distance helioseismol- 
ogy show the existence of large-scale motions of plasma in the 
sub-photospheric regions of the Sun (see, for example, a general 
review of recent findings in theoretical and observational loc al 
and time-distance helioseismology by iGizon & Birchl (|2(305)). 
These flows can be caused, for example, by interaction of 
the magnetic field with the plasma in the sunspots, by dif- 
ferential rotation, or by meridional circulation. An initial, en- 
couraging analyses of sub-photospheric flows in the limit of 
sho rt wavelength (ray approximation) has been carried out 
by iGilesI d!999l) . Time-distance helioseismology in the ray 
approximation has also been used to m a p the flow struc- 
beneath sunspots dKosovichev et al.L l2000t IZhao et ahl 



tures 

l200lt IZhao & KosovTchevl l2003h . The presence of long-lived 
large-sc ale subphot o spheri c conv ective flows has been stud- 
ied by lHaber et all (fl998l l200l using ring-diagram analy- 
sis. Surface gravity waves h ave been used in the studies of 
Corbard & Thompson (2002) to analyze radial gradients of an- 
gular velocity in the solar subsurface regions down to 15 Mm. 
Progress in the analytical analysis of solar sub-photosphere 
with a homogeneous ba ckg round flow has been shown by 
Erde lyi & Tarovanl d2001l) and lTarovanl d2004l) . 



The measurement accuracy of solar observations is contin- 
ually changing with the improvement of observational technol- 
ogy, thus requiring more precise methods of data analysis and 
giving a need for artificial (synthetic) data produced by mod- 
elling simulations, which are able to mimic observations with a 
required precision. The approach of forward modelling of sound 
speed inhomogeneities in homogeneous solar interior models 
was used to explor e the validity conditions of th e Born and ray 
approximations by Birch & Kosovichev (2000). Also, a study 
of the validity of the ray and Born approximations in analyz- 
ing the sensitivity of wave travel times to background flows, 
based on the forward modelling of acoustic wave propagation 
in a constant ba c kgrou n d density and pressure, has been done by 
iBirch & Felderi d2004t) . IGizon & Birchl d2005l) . Later, a limited 
approach using numerical solution of the wave equation with so- 
lar sound speed and density profiles has been applied to analyse 
seis mic traces produc ed by propagation o f acoustic waves, see 
e.g. lTong etail d2003l) . In a recent work of lShelvag elatl (2006) 
the ability of full forward modelling to successfully reproduce 
the observed power spectrum of solar oscillations and seismic 
traces by excitation of sound waves by single and multiple ran- 
dom acoustic sources has been shown. 

In this paper, an analysis is presented of the results of 
the forward modelling of acoustic wave propagation in the so- 
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lar sub-photosphere with two types of sub-photospheric steady 
horizontal flows. The full set of compressible hydrodynamic 
equations is solved in order to provide the artificial data. 
The initial model is calculated fr om the standard Model S 
dChristensen-Dalsgaard et all Il996l) . and is modified in order 
to suppress convective instability in the originally convectively 
unstable solar sub-photosphere. Flows with Gaussian and lin- 
ear velocity dependences on depth are analysed in the terms of 
the techniques of local time-distance helioseismology. The travel 
times and travel-time differences are computed from the time- 
dependent simulations and are compared with results calculated 
using the ray approximation. 

In the next section we give a brief description of the simula- 
tion setup and outline briefly the numerical methods used to pro- 
duce the artificial data. In Section 3, the simulations, their direct 
interpretation and data measurements in terms of time-distance 
helioseismology are shown. Further, the method of obtaining 
travel-time differences from the artificial data is discussed in this 
section. In Section 4, we compare the travel-time differences ob- 
tained from forward modelling with those calculated using the 
ray approximation. The inversion technique, the flow profiles ob- 
tained by inversion of the travel-time differences, and their com- 
parison with the original sub-photospheric flow velocity profiles 
are presented in Section 5. The results and conclusions are dis- 
cussed in Section 6. 



2. Simulation setup 



1000 km 



The s imulation model box has been described by lShelyag et alJ 
(2006). For convenience, we briefly describe the setup here. 
The Versa t ile Ad vection Code (VAC), originally developed by 
iTofh et al.l (Il998l) . is implemented here to carry out the forward 
modelling. VAC solves numerically a set of hyperbolic equations 
by a number of different computational methods in one, two or 
three dimensions in Cartesian, cylindrical or spherical geome- 
tries with uniform or non-uniform grids. The code was specif- 
ically set here to solve the ideal hydrodynamic equations in a 
two-dimensional Cartesian domain with three-dimensional vec- 
tor quantities: 



% + V • (vp) = 0, 
ot 



d(py) 
dt 



+ V ■ (vpv) + Vp = pg, 



de 

— + V • (ve + vp) = pg • v + Q (r, t) , 

Ot 



p = (IW) \e~ 



(1) 
(2) 
(3) 
(4) 



where p is the density, v = (v t , 0, v z ) is the two-dimensional ve- 
locity vector, e - e + pv 2 /2 is the total energy density per unit 
volume, e - pi (T\ - 1), p is the kinetic gas pressure, T\ is the 
adiabatic index, g = (0, 0, g) is the solar gravitational accelera- 
tion vector, Q is the temporally and spatially dependent term 
describing additional energy sources, and t is time. 

A Total Variation Diminishing (TVD) spatial and a fourth- 
order Runge-Kutta time discretisation schemes are used. The 
simulation domain is shown in Fig. 1. The box is 150 Mm wide 
and 50 Mm deep, and has a resolution of 960x4000 grid points; 
the upper boundary of the domain is near the solar temperature 
minimum (see Fig. [TJi. The boundaries of the domain are open 
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Fig. 1. The simulation domain. The measurement level is marked 
by the dashed line. The temperature profile is also sketched on 
the plot. 
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Fig. 2. Dependence of the gradient of the sound speed on depth 
in the domain. 



with zero gradients across the boundaries. Two boundary regions 
with g = at the top and bottom boundaries of the box are intro- 
duced in order to simultaneously keep the domain in the pressure 
equilibrium and to satisfy the boundary conditions. 

The equation of state in the code is used in the form, which 
describes the properties of an ideal gas. The adiabatic index 
Fi, assumed to be constant over the computational domain, is 
equal to 5/3. The pressure equilibrium and modified convective 
stability conditions with constant solar g ravity accelera t ion ar e 
applied to calculate the other quantities (Shely ag et all I2006). 
The upper density value and convective instability dependence 
(Schwartzschild criterion) taken from Christensen-Dalsgaard's 
standard Model S, are used as the initial conditions for the pres- 
sure equilibrium calculation. A slight increase of F; in compar- 
ison to the original solar model makes most of the domain ex- 
cept the super-adiabatic strongly convective layer near the solar 
surface convectively stable. Thus, the convective instability de- 
pendence is modified in this layer in order to reach convective 
stability. This procedure of pressure equilibrium recalculation 
results in values of sound speed close to those in the Sun. The 
dependence of the sound speed gradient on depth in the domain 
is shown in Fig. [2] 

The perturbation source is located in the upper-middle 
(z=2000 km below the upper boundary) of the simulation box 
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Fig. 3. Time-distance diagram of the acoustic response of the 
model to the source. The first and second bounces are marked 
approximately by dashed and dash-dotted lines, respectively, and 
are clearly visible. 

(see Fig.[T|i. The measurement level for latter time-distance anal- 
ysis is located at 1000 km below the upper boundary. 
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Fig. 4. Sub-photospheric flow structures introduced to the sim- 
ulation domain. Dashed line corresponds to the linear flow ve- 
locity profile with the maximum speed of 100 m s _I at 50 Mm 
depth and zero at the source location, solid line corresponds to 
the Gaussian flow profile with the maximum of about 330 m s _1 
and the FWHM of 10 Mm. 



3. Simulations 

In order to excite so und waves, a single perturbation source (see 
Shel yag et all [2006). corresponding to a localised cooling even t 
causing mass inflow and sound waves excitation dRasti [1999), 
has been introduced in the simulations. The source, dQ, is de- 
scribed in the energy equation of system of hydrodynamic equa- 
tions as 



dQ = 



1 + tanh log 3 



t-tp 



■ exp ■ 



( r-rp) 1 

CT2 2 



(5) 



The characteristic timescale of the source is cr\ — 120 s, 
the spatial extent <T2 = 0.5 Mm. The source is located in the 
middle of the horizontal layer at 2000 km depth of the simulation 
domain. 

The acoustic response of the model is shown in Fig. [3] The 
plot shows time-distance diagram of the vertical velocity com- 
ponent taken at the measurement level of 1000 km below the 
upper boundary of the simulation domain. This level approxi- 
mately corresponds to the level of the visible solar photosphere. 
The first (faster) and the second (slower) bounces of the sound 
wave, propagating in the simulated solar interior, are visible in 
the plot. Also, the slowest wave, propagating in the waveguide 
created by the simulated temperature (and sound speed) mini- 
mum, is seen in the figure. 

We study the influence of steady flows to the acoustic 
response of the simulated solar sub-photosphere. Two types 
of flows are introduced (see Fig. |4|. The first one is a lo- 
calised flow with Gaussian horizontal velocity profile. The 
flow centre is located at the depth of 25 Mm. The width 
of the flow (FWHM) is 10 Mm. The maximal flow veloc- 
ity is about 300 ms _1 . The flow can represent, for exam- 
ple, horizontal motion of solar plasma near a large sunspot. 
Existence of such mass flows beneath sunspots has been dis- 
covered by the means of time-distance helioseismology by 
Zhao et al. 1 (120011) using the ray approximation. In this work it 
has been shown that localised horizontal subsonic flows with the 



speeds of the order of kilometre per second exist at depths of 
about 10 Mm and have the horizontal extent of about 30 Mm 
dZhao et al 1 120011 Fig. 3 in their paper), the results of inversions 
suggest approximate cylindrical symmetry of the flow struc- 
tures under sunspots. Full magneto-hydrodynamic simulations 
of magneto-convection in convectively unstable cylindrically- 
symmetric polytropi c model with presence of s trong magnetic 
flux concentrations (Hurlbur t & Ruckli dge. 2000) also suggest 
the existence of such sub-photospheric large-scale convective 
flows. These flows in the simulations appear essentially weakly 
magnetised in comparison to strong magnetic flux concentra- 
ti ons, where convection is sup pressed by magnetic field (Fig. 2 
in lHurlburt & Rucklidgd 12000) . 

The second case is characterised by a linear velocity depen- 
dence on the vertical coordinate selected in the way that the ve- 
locity at the source location is equal to zero, and the velocity 
at the bottom of the domain is 100 m s _1 . This flow can model 
meridional lam inar circulation. Such fl ows have been seen obser- 
vationally bv lHaber et al.l dl998ll20Qlb using ring-diagram anal- 
ysis of SOI-MDI observational data. Simulations of convection 
in spherical shells also show an evidence fo r meridional circu- 
lation flows with su fficiently long lifetimes (Ell iott et all 12000; 
iMiesch et al.1 12000). The average dependences of the flow speed 
on depth for different latitudes and C arrington longitude s, ob- 
tained by inversions of ring diagrams (lHaberetall lT998. their 
Figs. 5, 6), have qualitatively similar properties as the linear 
flow introduced to the numerical experiments presented here. 
The flow profiles are approximately linearly dependent on depth, 
and the speeds at the depth of 15 Mm are of 25 - 50 m s . 

Since the acoustic source is equidistant from the side bound- 
aries of the simulation box, and without a flow the wave packets 
would propagate symmetrically around the source, one can study 
effects, caused by the Doppler shift differences between the wave 
packets propagating to the left and to the right from the source. 
The numerical difference between the left and right amplitudes 
represents the wave phase shift. 

The phase difference image for the case of the Gaussian flow 
profile is shown in Fig. [5] which shows the difference between 
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Fig. 5. Vertical speed difference image for the model with the 
Gaussian horizontal velocity profile. The differences are com- 
puted between the points located at the same distance and oppo- 
site sides of the source. Positions of the first and second bounces 
in the diagram are marked by dashed and dash-dotted lines, re- 
spectively (cf. Fig. 3J. The figure shows that only the first (the 
deepest propagating) bounce is influenced by the flow. Also, the 
flow acts on the wave propagating in the deeper regions with al- 
most constant sound speed (Fig. [2J, thus, the "difference waves" 
travel with a constant speed, corresponding to the sound speed in 
the region of the model, where the influence of the flow is most 
significant. 

the right and left portions of a time-distance diagram, such as 
Fig. [3] The introduced localised flow acts only on the modes 
which propagate deep enough in the solar sub-photosphere; the 
shallow modes are much less affected. This fact can be clearly 
seen from the figure. The phase differences are large in the re- 
gions of the first bounce propagation, however, the second and 
higher order bounces and the surface wave are not influenced by 
the flow. Also, since the flow is localised in the deep regions of 
sub-photosphere, where the sound speed is changing slowly, the 
phase difference "waves" are propagating as the waves with con- 
stant group speed corresponding to the local sound speed at the 
region where the influence of the flow is most significant, as can 
be seen from the figure. 

The case for the linear horizontal speed profile is rather dif- 
ferent. The phase difference plot calculated for the linear profile 
is shown in Fig. [6] In this case all of the modes are influenced by 
the flow. 

4. Travel-time differences 

We perform the comparison of the simulated travel times with 
the travel times given by theoretical calculations using the ray 
approximation for the model depth dependences of the sound 
speed, pressure and density used in the numerical experiments. 
The simulated travel-time differences are calculated from the 
phase differences between the wave packets propagating to the 
left and to the right from the source. Assuming the harmonic re- 
sponse of the model to the source, the signals to the left and to the 
right of the source can be defined as A/ = Aq sin (at (t - 6t/2)) 
and A r - Aq sin (lo (t + St/2)), where Aq is the amplitude of 
the wave packet, t is time, and St/2 is the phase shift of the 
wave packet with respect to the wave packet propagating in 
the model, undisturbed by the flow. Thus, assuming that the 
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Fig. 6. Vertical speed difference image for the model with the 
linear horizontal velocity profile Fig. [5] All modes in the box are 
influenced by the flow. 



phase differences are less than 2n and measuring the difference 
AA - A\ - A r - 2Ao cos (ut) sin {w ■ St 12) from the simulations 
accordingly to Figs.|5]and|6]one can evaluate the time difference 
between the acoustic waves propagating to the left and to the 
right from the source: 



2 AA/A 
St — — arcsin — - — . 
co 2 



(6) 



In Eq. ©, Aq is the amplitude of the oscillation with the fre- 
quency a), taken along the time-distance path of the wave packet, 
AA is the amplitude difference between the packets propagating 
to the left and to the right from the source, taken at the same 
place where the amplitude Aq is taken. 

In order to calculate the travel-time differences i n the ray ap - 
proximation, we follow the procedure described bv lGilesI |l999), 
appropriately adapted the method for the Cartesian geometry 
used in the simulations. At the solar surface, travel-time differ- 
ence between the wave packets propagating to the left and to the 
right from the source can be written in the form: 



St 



K ■ u dz, 



(7) 



where St is the travel-time difference, u is the horizontal flow ve- 
locity, dependent on the depth, K is the integration kernel, which 
contains the information about the physics of the wave propaga- 
tion depending on the depth. The integration limits (zi , Zi) are 
defined by the depth range where an oscillatory mode of some 
frequency a> can propagate, or, in other words, by the lower (zi) 
and upper {zi) turning points for this mode. 

In the ray approximation theory, the kernel K is defined by 
the following expression: 



K = 4- 



(8) 



The factor 4 in Eq. (O indicates that in the case of the simu- 
lations described above, the travel-time differences for the waves 
propagating to the left and to the right are the same, and that the 
time the acoustic mode travels from the upper turning point to 
the lower one is equal to the travel time in the opposite direc- 
tion. 
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In a plain-parallel model one can write the quantities under 
the integral as follows: 



k 7 a> c 



3 „2 



V 7 = 



cj a - k 2 c 2 cj 2 BV ' 



krLOC 



2 2 
2 w -°>BV 



■ k 2 C 2 CJ 2 



where cobv is Brunt-Vaisala frequency, defined as 
2 /Winn dlnp\ 

o»ac is the acoustic cutoff frequency, given by 



c l / din p 



dz 



k z is vertical wavenumber, 



c 2 ( a) 



mac) 



(9) 



(10) 



(11) 



(12) 



(13) 



and, finally, k x is horizontal wavenumber. 

The upper and lower turning points zi and Z2 are defined as 
the range where the integrand of Eq. (|7]) is real. 

Substituting Eqs. (|9]l-(fT3l into Eq. ([8]), one can calculate the 
values of the kernel, defined by Eq. (|S). However, numerical 
computations of the integral in Eq. (0 with this kernel around 
the lower turning point has to be carried out carefully due to 
the singularity of the integrand there. This singularity is inte- 
grable; the procedure of integratio n is described in details by 
IChristensen-Dalsgaard et al. I ( Tl989l) . 

The results of the comparison are shown in Figs.|7]and[8]for 
the Gaussian and linear horizontal velocity profiles, respectively. 
In the case of the Gaussian velocity profile (Fig. [TJ, the ray ap- 
proximation does not give correct travel times for medium (30- 
60 Mm) distances. However, the difference between travel times 
derived between the ray approximation and the simulations di- 
minishes on larger distances for the waves that propagate deep 
in the sub-photosphere. 

The dependence of the travel-time difference on the horizon- 
tal coordinate for the linear velocity profile shows a rather dif- 
ferent behaviour. All of the sound waves propagating in the sub- 
photosphere are affected by the flow, and the difference between 
the ray approximation and the forward simulations grows with 
the distance from the centre. Also, one can notice that the trend 
of the dependences show qualitatively similar behaviour. 

The noise in the travel-time difference dependences at dis- 
tances about 60-70 Mm from the source, is caused by the inter- 
ference of the forward wave packet with the wave reflected from 
the boundary of the simulation domain, despite the reflected 
wave, having much lower amplitude than the one propagating 
toward the boundary. 

Generally, it seems that the ray approximation does not give 
accurate enough travel times when compared with the ones de- 
rived from the simulations for the same sub-photospheric struc- 
ture. This is caused by the fact that the wave packets, produced 
by superposition of p modes, propagate along the ray bundle 
which follows the path given by the ray approximation, but has 
a finite extent (Bogdan, 1997). Thus, for the case of a Gaussian 
flow, the wave packets actually reach deeper layers of the sub- 
photosphere and become more affected by horizontal flow, then 
predicted by the ray approximation. This difference diminishes 
when the turning point of the ray path is located close to or 
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Fig. 7. Plot of travel-time difference for "left" and "right" prop- 
agating wave packets in the case of localized Gaussian flow pro- 
file. The solid curve corresponds to the time difference measured 
directly from the simulations, the dash-dotted curve is the time 
difference calculated using the ray approximation for density and 
pressure profiles used in the simulations. The ray approximation 
does not seem to give accurate results for small distances (and 
depths) in this case, however, for large distances the difference 
diminishes. The noise at distances about 70 Mm is caused by in- 
terference of the forward-propagating wave packet with the wave 
reflected from the boundary of the simulation box. 
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Fig. 8. Same as Fig. [7] computed for the model with linear ve- 
locity depth dependence. In this case the flow is not spatially 
localised, and the time difference grows with the distance from 
the source. 



below the maximum of the horizontal flow. The wave packets, 
travelling in the sub-photosphere with a linear horizontal flow 
with flow speeds increasing toward the bottom of the simulation 
box, again because of their finite extent, are more influenced by 
the deeper and stronger flows in comparison with those of the 
ray approximation. However, as it has already been noticed, the 
travel-time difference dependences are qualitatively similar: for 
example, in Fig. [8]at e.g. the distance x around 10-20 Mm, both 
dependences show a similar feature. This is a change of incli- 
nation of the curves, which is connected to the change of sound 
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speed gradient (cf. Fig.[2| around the depth 5-10 Mm in the sub- 
photosphere. This change appears at smaller distances for the 
simulations in comparison with the ray approximation, because 
the real wave packets reach the region of sub-photosphere with 
a smaller sound speed gradient earlier than the wave packets ap- 
proximated by ray theory, due to their finite extent. 

5. Inversions of velocity profiles 

Next, we perform inversion of the flow velocity profiles in the 
domain using the ray approximation. The i nversions are carried 
out following the procedure described by Giles] d!999t) . In the 
case of spatially discretized grid, Eq. (0 can be cast in a matrix 
form: 



K ■ u — 6t, 



(14) 



where r = r, is the vector of travel-time difference measure- 
ments, u - uj is the vector representing the flow speed depen- 
dence on depth, K is the matrix representation of the integration 
kernel Eq. (|S): 



(15) 



The integral is taken along that part of the ray path which lies 
within the grid element ij. 

The inversion problem can be set as a proble m of rmnirnisa - 
tion of the quantity \Ku - £t| 2 (see, for example. iPressl 12002). 
which is equivalent to 



K T Ku = K t 6t. 



(16) 



In general, the problem described by Eq. ( fTSI l can be ill- 
posed, also the setup can deal with the noisy data, thus the opera- 
tor K T K has to be regularised. In order to diminish the numerical 
noise coming from the simulations, we regularise Eq.(fT6b in the 
following way: 



(k t K + yh) u = K t 5t, 



(17) 



where A is the regularisation operator, and y is a parameter 
which controls the relative influence of regularisation on the 
minimisation procedure. 

For the inversion problem, the system of equations Eq. (JT7J 
is solved with respect to the flow velocit y profile u using the 
singular value decomposition dPressl 120021) . 

The forward problem of calculating 6t, based on the given 
velocity profile and model structure, has been solved using the 
ray approximation. Then, the obtained values of 6t have been 
used in the inversion process in order to test the method, to select 
the regularisation parameter y, and to compare the flow profiles 
obtained from the simulations and from this calculation, consis- 
tent with the ray approximation. The regularisation parameter 
y in Eq. (jT7j is selected in such a way that it does not signifi- 
cantly change the velocity profile, obtained by inversion of the 
travel-time difference dependence, calculated using the ray ap- 
proximation, in comparison with the original one. The linear reg- 
ularisation operator of the form, describing a priori solution as a 
constant, is chosen for the inversion procedure. 

The results of inversion are presented in Figs. [9] and [10] Solid 
curves in these plots correspond to the original flow profiles, 
the dashed ones are obtained from the inversions of the travel 
time differences, which have been calculated using the ray ap- 
proximation for the original flow profiles, and the dash-dotted 
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Fig. 9. Inversion of the velocity profile from the travel-time dif- 
ference measurements using the ray approximation. The solid 
line corresponds to the original flow profile. The dashed line is 
the velocity profile calculated from the travel-time differences 
obtained by integration along the ray path given by the ray ap- 
proximation theory for the model thermodynamic parameters. 
The dash-dotted line shows the velocity profile inferred from the 
travel-time difference measurements of the simulations using the 
ray approximation. The profile, inferred from the inversion, ap- 
pears to be broader and has lower maximal amplitude than the 
original one. 
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Fig. 10. Same as Fig. [9] calculated for the model with linear ve- 
locity dependence on depth. The inferred profile shows larger 
flow velocity values than the original one at depths above ap- 
proximately 25 Mm, however, it starts to decrease in the deeper 
sub-photospheric layers because of the geometry of the model. 



curves show the flow profiles which are obtained from the in- 
versions of the simulated travel-time differences (Figs. [7] and [8] 
solid curves). 

In the region with approximate depth of less than 30 Mm 
from the solar surface, the original profile (Figs.l9land[T0l solid 
curves) and the profile obtained by inversion of the travel-time 
differences, calculated from the ray approximation (Figs. [9] and 
[TOl dashed curves), correspond well to each other. This defines 
the validity range of the data inferred from the inversion proce- 
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dure. In the both flow cases, the results for the regions with the 
depth more than about 30 Mm do not represent the flow pro- 
files precisely, due to the geometry of the numerical setup. In 
the model setup the wave packets that reach the deeper regions 
come out of the domain through the side boundaries and do not 
reach the surface. 

The velocity profiles, obtained by inversions of the travel- 
time differences, calculated from the simulations (Figs. [9] and 
[TUl dash-dotted curves), show significantly different behaviour 
compared to the real flow profiles. In the simulations, the wave 
packet, while propagating in the sub-photosphere, is influenced 
by the deeper regions of it in comparison with the ray approx- 
imation. Thus, in the case of horizontal flow with the Gaussian 
depth dependence, the profile, obtained from the inversion, ap- 
pears to be broader, and has lower maximal amplitude than 
the original one, because the sensitivity of the wave packet is 
smoothed over the range of depths. Due to a similar reason, the 
velocities, obtained by the inversion, in the region between 20 
and 10 Mm are larger, than the velocities in the original profile. 

In the case of the linear flow velocity profile, the velocity 
dependence, obtained by inversion, behaves in a rather differ- 
ent way. Above approximately 25 Mm depth level, the obtained 
profile shows velocities larger than the velocities of the original 
flow. Below 25 Mm the dependence starts to decrease toward the 
bottom of the model. Again, this is caused by the finite width of 
the wave packet. A part of the wave packet leaves the domain 
after it reaches the lower boundary. 

Generally, the both flow profiles are consistent with the cor- 
responding travel-time difference dependencies Figs. [7] and [8] 
Larger travel-time difference corresponds to larger influence of 
acoustic wave by the flow, which means the larger flow speed at 
the same depth level. For the deeper regions, the decrease of the 
flow speed with depth is partly due to the geometry of the simu- 
lation setup, and in the case of Gaussian flow profile the decrease 
is also caused by the localised nature of the introduced flow. 

The width of the wave packet at a particular depth can 
roughly be estimated from Fig. [10] The width should not be less 
than the distance difference between the points on the curves of 
the profile obtained from inversion of the simulated data and the 
profile obtained using the ray approximation taken for the same 
value of the velocity of the flow. 

6. Discussion 

In this paper, analysis of the influence of sub-photospheric flows 
on acoustic wave propagation in the solar interior is presented. 
The numerical code implemented here solves the full set of com- 
pressible hydrodynamic equations in two dimensions. The ini- 
tial model uses constant adiabatic index Fi = 5/3 and solar 
Standard Model S modified in order to suppress convective in- 
stability in the solar interior, which results in values of the sound 
speed close to those in the Sun. Acoustic waves are generated 
by a source corresponding to photospheric cooling events (pho- 
tospheric plumes). 

Using forward modelling, the travel times and travel-time 
differences have been obtained and compared with the travel- 
time differences calculated using the ray approximation tech- 
nique for two different sub-photospheric flow profiles: localised 
horizontal flow with Gaussian dependence of the flow speed 
on depth, and flow with a linear dependence of the flow speed 
on depth. The Gaussian flow may represent sub-photospheric 
plasma motion surrounding a sunspot. The linear flow may 
mimic large-scale flows. The comparison of the travel-time dif- 
ferences shows that, in the case of a Gaussian equilibrium mo- 



tion centered about 20-30 Mm below the solar surface, the travel- 
time difference between the ray approximation and forward 
modelling is large at small distances from the acoustic source 
and diminishes at large distances. The behaviour of acoustic 
waves propagating in a model with linear velocity profile is dif- 
ferent: there the travel-time difference grows with the distance 
from the source. This difference can be explained by including 
some wave effects, which are neglected in the ray approxima- 
tion. Also, the discrepancy between actual travel times and those 
given by the ray approximation creates certain limitations which 
have to be taken into account in analysis of real observational 
data. 

The influence of the wave effects has been analysed by inver- 
sion of the travel-time differences obtained for both cases of the 
sub-photospheric flow velocity profiles. The inversion of travel- 
time differences was carried out using the ray approximation. 
Then the dependences of the horizontal flow velocity on depth, 
evaluated by inversion, were compared with the original equi- 
librium bulk motion profiles. The comparison shows that in the 
depth range between and 20 Mm, the velocities obtained by 
inversion of the modelled data appear to be larger than the orig- 
inal ones; however, in the case of Gaussian flow, the maximal 
velocity is lower for the profile obtained from inversions. This is 
caused by the sensitivity of the wave packet, which in reality is 
spread over a range of distances, while in the ray approximation 
the wave packet is presumed to be infinitely thin. 
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